Template coexistence in prebiotic vesicle models 
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The coexistence of distinct templates is a common feature of the diverse proposals advanced to 
resolve the information crisis of prebiotic evolution. However, achieving robust template coexistence 
turned out to be such a difficult demand that only a class of models, the so-called package models, 
seems to have met it so far. Here we apply Wright's Island formulation of group selection to study 
the conditions for the coexistence of two distinct template types confined in packages (vesicles) of 
finite capacity. In particular, we show how selection acting at the level of the vesicles can neutralize 
the pressures towards the fixation of any one of the template types (random drift) and of the type 
with higher replication rate (deterministic competition). We give emphasis to the role of the distinct 
generation times of templates and vesicles as yet another obstacle to coexistence. 



I. INTRODUCTION 



Though a ubiquitous feature of life (e.g., genomes, 
multi-cellular organisms and animal societies) coopera- 
tion poses a difficult problem to the classical interpreta- 
tion of evolution by natural selection as an egoistic pro- 
cess in which individual organisms compete fiercely to 
guarantee the presence of their genes in future genera- 
tions. Cooperative traits are costly because they presup- 
pose the investment of resources towards a public good, 
thus benefiting other individuals who may utterly fail to 
contribute to the community welfare. Traits that benefit 
the group as a whole but are deleterious to their bearers 
used to be termed altruistic but, probably to avoid the 
heavy anthropocentric connotation of this term, the mod- 
ern literature favors the denomination cooperative traits 
instead There is an extensive literature on the evo- 
lution of the cooperation in nature (see e.g. [H, 0> H, 0| ) 
but only very recently a series of experiments on micro- 
bial populations have substantiated the theoretical argu- 
ments supporting it @, H 0| • 

A key element to explain the evolution of cooperative 
traits is the existence of some form of non-random as- 
sociation between individual members of the population. 
For instance, cooperation may be restricted to a group 
of individuals sharing a common ancestor. This is the 
so-called kin selection mechanism that gained fame by 
explaining the cooperative behavior of social insects [8(. 
Curiously, although kin selection is fittingly considered 
a kind of group selection, it relies entirely on a gene's- 
eye view in which the unit of selection is the gene rather 
than the individual or the group. A robust alternative 
to the postulation of mechanisms of recognition between 
kin individuals is the spatial enclosure or compartmen- 
talization of small groups of individuals since then the 
offspring of cooperators will remain close to their rela- 
tives and the benefits from cooperation will be confined 
mainly to the group of cooperators. In fact, compartmen- 
talization is now acknowledged as an obligatory stage in 
prebiotic evolution needed to enforce cooperation among 
distinct templates and so to set the conditions for the 
formation of a gene network necessary to cellular life @ . 



Henceforth, we will use interchangeably the terms tem- 
plate and replicator to refer to a self-replicating molecule 
formed by a sequence of nucleotides. 

The awareness of the necessity of a primitive mecha- 
nism to impose cooperation among templates has grown 
from the seminal work of Eigen on purely competitive 
templates which came to be known as the quasispecies 
model 10]. Particularly relevant to our purposes was the 
finding that the length of a replicating polymer (i.e., a 
RNA-like template) is limited by the replication accuracy 
per nucleotide and so primordial replicators would have 
to replicate with very high accuracy in order to reach 
the length of today's RNA viruses (about 10 3 to I0 4 nu- 
cleotides). Such a replication accuracy, however, cannot 
be achieved without the aid of specialized catalysts (pep- 
tide enzymes), but building those catalysts requires a 
blueprint that amounts to a large nucleotide sequence, 
leading thus to a molecular version of the old puzzle 
about the chicken or the egg - no large genome without 
enzymes, and no enzymes without a large genome fTTj | . 
This finding and the observation that templates with dis- 
tinct replication rates cannot coexist [12( triggered the 
so-called information crisis of prebiotic evolution. This 
crisis brought to light the challenging conclusion that, in 
spite of being at the heart of natural selection, compe- 
tition alone would not have worked in prebiotic times: 
some form of cooperation between templates is manda- 
tory to overcome the information crisis. 

Cooperation can bypass the information crisis by al- 
lowing the coexistence of several short templates, i.e., 
by splitting the information in modules, similarly to 
the division of the genome in chromosomes found in 
many organisms. To ensure the coexistence of distinct 
replicators, Eigen and Schuster proposed a cyclic reac- 
tion scheme, termed hypercycle, in which each replicator 
would aid in the replication of the next one, in a regu- 
latory cycle closing on itself [13]. Clearly, this proposal 
requires that the primordial replicators functioned both 
as template and replicase. Although the discovery of the 
catalytic activity of RNA has lent credibility to the hy- 
percycle scenario [T3], the assumption that each replica- 
tor has two separate functions, namely, a replicase for 
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the next member of the hypercycle and a target for the 
replicase associated to the previous member faced strong 
criticism [jq . In fact, while it is obvious that natural 
selection will act so as to make each element of the hy- 
percycle a better target for replication, it will oppose or 
at least not favor the cooperative part of the scheme, i.e., 
to make the replicator a better replicase for other replica- 
tors. This function is then certain to degenerate quickly 
since deletions and mutations that impair it would carry 
a selective advantage. In that sense, as Maynard Smith 
pointed out, giving catalytic support in such molecular 
networks is an altruistic behavior and so hypercycles are 
easy targets to parasites, i.e., molecules that do not recip- 
rocate the catalytic support they receive (l6| . The ruin 
of the hypercycle is then an unavoidable consequence of 
natural selection. 

An alternative suggestion to resolve the problem of the 
coexistence between templates which is very much in line 
with the classical works on the origin of life [l?], [H| is to 
enclose the unlinked templates in isolated compartments 
or vesicles. The key to coexistence is to assume that 
the vesicles proliferate with a production rate that de- 
pends on their template compositions. Essentially, this 
amounts to assume that the coupling among different 
template types occurs through a common metabolism 
which is ultimately the responsible for the survival and 
reproduction of the vesicle, and that the well-functioning 
of this metabolism requires the contribution of all tem- 
plate types. This is the central idea behind the so-called 
package model s fl9l. a mong which the stochastic correc- 
tor model (20I . l2lL I22L I23I ] is the most popular. As re- 
vealed by the word 'stochastic', a crucial ingredient of 
these models is the finitude of the population of tem- 
plates within each vesicle. In fact, it is the stochastic 
nature of the template dynamics that produces diversity 
in the population of vesicles, creating thus the opportu- 
nity for the operation of natural selection. 

In contrast with the quasispecies and hypercycle mod- 
els, very little is known about the dynamics and station- 
ary states of package models, since the great complex- 
ity resulting from the coupling of template and vesicle 
dynamics precludes a full analysis of the space of pa- 
rameters of the models. Such a systematic analysis is 
important to determine in what conditions, if any, tem- 
plate coexistence can be achieved. The situation here is 
similar to that found in models of group selection (see, 
e -g-, [1, H> 0] for reviews). In particular, the coupling 
between the two dynamics can be treated analytically 
provided there is a countable infinity of vesicles, so the 
dynamics at the group level is deterministic [U, [25|, H|| . 
Otherwise, this group selection formulation retains the 
main ingredients of the package models in that the num- 
ber of templates within each vesicle is finite and the sur- 
vival and consequent proliferation of the vesicles depends 
on their template compositions. In this contribution we 
broaden a preliminary study on the suitability of the 
classic group selection framework to study coexistence 
of templates in package models [27j . In particular, we re- 



lax the unfounded but widely used assumption of group 
selection models (see, e.g., [2J, [2|| [2|| H3] ) that templates 
and vesicles have the same generation times. In doing so 
we found that coexistence is impossible if the template 
generation time is much shorter than that of the vesicles 
and that, when possible, template coexistence is achieved 
only within a well-defined range of the vesicle capacity. 



II. MODEL 

Following the formulation of Wright's Island model [28[ 
we consider a global population composed of an infinite 
number of spatially isolated local populations - the vesi- 
cles - each of which encloses exactly N templates. This 
framework, which forms the foundation for traditional 
group-selection theories, has been successfully used to 
study the conditions for the evolution and maintenance 
of altruistic traits in nature (see, e.g., (H). Here we em- 
ploy the Island model to study a more difficult problem, 
namely, the coexistence of two distinct template types A 
and B within a same vesicle of capacity N . Without loss 
of generality we assume that template A has replication 
rate 1 — r with t £ [0, 1] and template B replication rate 
1. Hence r is referred to as the handicap parameter. The 
vesicles are identified by their template compositions or, 
more pointedly, by the number i = 0, . . . , N of type A 
templates they wall in. The state of the global population 
at a given generation t is completely specified by the fre- 
quencies of vesicles of type i (i.e., a vesicle that encloses i 
templates A and N — i templates B), denoted by Yf with 
Y* — 1 for all t. Given the mechanisms for template 
competition that takes place inside the vesicles and for 
the competition between vesicles, our goal is to derive a 
recursion equation for the frequencies of the N + 1 differ- 
ent vesicle types. The life cycle of the vesicles (i.e., one 
generation) consists of three events - vesicle extinction, 
vesicle recolonization and template replication - which 
take place in this order and are described as follows. 



A. Vesicle extinction 

As pointed out in Sect. HI the coexistence of distinct 
template types solves the information crisis problem be- 
cause the information content of a vesicle may be seen 
as split into several (two, in our case) parts, which must 
all be present at any time to guarantee the viability of 
the system as an integrator of information. Moreover, it 
would be desirable that the different templates contribute 
with approximately the same number to the vesicle com- 
position. This constraint, namely, the presence of both 
template types in equal concentrations within each vesi- 
cle, can be enforced by choosing an appropriate measure 
for the survival probability of type i vesicle, which we de- 
note by cti £ [0,1]. Here we choose the geometric mean 
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(1) 



where g G [0, 1] is a parameter measuring the benefit 
brought to the vesicle by the presence of the two tem- 
plate types within it. Hence regardless of the value of 
g > survival is guaranteed (ai = 1) for vesicles with an 
even template composition (i — N/2). Vesicles that lack 
one of the template types, i.e., vesicles of type i = and 
i = N are assigned a baseline survival probability 1 — g, so 
that the selective pressure against these vesicles increases 
with increasing g. We note that the vesicle survival prob- 
ability given in equation |T]) can be seen as describing 
a dynamical link among the template types through a 
common metabolism - each template contribute to the 
good of the vesicle by catalyzing its metabolism at var- 
ious points [2(J. The absence of any of the catalysts 
would then greatly impair the vesicle metabolism and so 
its survival capability. 



B. Vesicle recolonization 

The net result of the extinction procedure described 
before is that a fraction 1 — ^ - a^Yi of vesicles disappear 
and must then be recolonized, i.e., replaced by the sur- 
viving vesicles. This is done by replicating these vesicles 
in proportion to their frequencies in the population just 
after the extinction procedure, yielding thus the following 
new vesicle frequencies 



(2) 



for i = 0, . . . , N. This equation prompts the interpreta- 
tion of the parameter cti as the replication rate of a vesi- 
cle of type i and so henceforth we will refer to the joint 
processes of extinction and recolonization as the process 
of vesicle replication. We note that this standard pro- 
cedure for recolonization (see, e.g., [IS US US] ) implies 
the instantaneous replacement of the extinct vesicles by 
the surviving ones with the probability given in equation 
(|2|). For example, in an extreme situation, in which only 
one vesicle passes the extinction stage, this sole vesicle 
will replenish the entire population (infinite or finite) in 
a single time step. Although rarely made explicit, this 
assumption seems to underlie all deterministic popula- 
tion models with non-overlapping generations. However, 
this drawback can be safely ignored if g is not close to 1 
which, fortunately, is the relevant situation in prebiotic 
evolution (see discussion in Sect. IVip . 



C. Template replication 

Since the capacity of the vesicles is fixed and finite, it 
is necessary to use a stochastic approach to model the 



dynamics of the templates inside each vesicle. As usual, 
we assume that the number of copies that a template 
brings forth is proportional to its relative replication rate. 
In addition, assuming that there is no overlap between 
consecutive generations of templates, i.e., between the 
parent and its clones, we can write the probability that 
a vesicle of type j changes to a vesicle of type i as Rij = 
(T™ 1 ) ^- where the elements of the matrix T are given by 
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(3) 



and wj — j(l — t)/(N — jr) is the relative replication 
rate of templates of type A in a vesicle of type j. Clearly, 
Y) . Rij = 1 Vj. The new ingredient here is the integer pa- 
rameter to > 1 that yields the number of replication cy- 
cles each local template population goes through for each 
generation of the vesicle population. In other words, if we 
set the generation time of the vesicles to 1, then the gen- 
eration time of the templates will be 1 /m, and so we refer 
to to as the ratio between vesicle and template genera- 
tion times. Our formulation for the template dynamics is 
essentially the celebrated Wright-Fisher model of popu- 
lation genetics [29| , which is very well suited to describe 
stochastic effects (e.g., random drift) in finite popula- 
tions. Up to now only the extreme case to = 1, in which 
templates and vesicles have the same generation time, 
was considered in the literature 0, [5S US H3| , though 
it is clear that a more plausible scenario would be to con- 
sider to 3> 1. We note that in some alternative prebiotic 
package models, in which the size of the local template 
population is allowed to increase, the vesicle replication 
stage is triggered when the template population reaches a 
certain size which is fixed a priori [IS, 0, [2^, US] . Thus, 
fixing this limiting size is similar to fixing our parameter 



D. Dynamics 

Finally, given the events comprising the life cycle of 
templates and vesicles we can write a recursion equation 
for the frequencies of vesicles of type i = 0,1, ... ,N in 
the global population 



yt+l _ 1^3=0 tMjWjIj 



E 



(4) 



Even in the stationary regime t — > oo, a closed solution 
for Y* is not possible, except when the template and 
vesicle dynamics decouple, which happens for g = and 
to — > oo. In both cases random drift ensures that for 
finite N either template A or template B will reach fix- 
ation within a given vesicle, i.e., 3^°° = for i ^ 0,JV. 
In fact, the very possibility of fixation of the less fit tem- 
plate A in a few vesicles together with the increment of 
the survival probability of those vesicles are key ingredi- 
ents of the classic models for the evolution and stability 
of altruistic traits [H, [H . 
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FIG. 1: Semi-logarithm plot of the stationary frequency 
of vesicles containing the less fit template as function of 
the handicap parameter for independent vesicles of capac- 
ity N = 6, 12 and 24. The solid lines are the exact results 
obtained from iteration of the recursion equations and the 
dashed lines are the results of the diffusion approximation. 



We note that our approach, based on Wright's frame- 
work for spatially structured populations, contrasts with 
Wilson's formulation Q in which the group structure is 
dissolved each generation to form a global mati ng p ool 
(see [30|, HH for application in ecology and [H, [33| for 
application in prebiotic evolution). In particular, coexis- 
tence is favored in this transient group formulation since 
different templates in distinct vesicles are likely to be as- 
signed to the same vesicle during the group re-assembling 
procedure after the mating stage. 

In this contribution we do not take into account the 
possibilities of interchange of templates between vesicles 
(migration) and mutations that change template A into 
B and vice-versa. These processes actually promote co- 
existence by preventing the fixation of the templates and 
so they are important to test the robustness of models for 
the evolution of altruism, the aim of which is the fixation 
of the less fit template A rather than the coexistence 
between the two template types [26|. In the following 
sections we will characterize the stationary solutions of 
the recursion equations (0]) for a wide range of the control 
parameters of the model. 



III. INDEPENDENT VESICLES 

In order to better appreciate the many obstacles hin- 
dering the coexistence of distinct templates confined in 
vesicles of finite capacity, in this section we focus on 
the simpler problem in which the vesicles evolve inde- 
pendently of each other. As pointed out before, this is 
achieved by setting g = in equation ([]} so that the 
survival of the vesicles is guaranteed regardless of their 
template compositions. Alternatively, by letting m — > oo 
we allow random drift to fix one of the templates (i = 
or i = N) and since ao = ajy = 1 — g the competition 
between vesicles is effectively turned off. Actually, the 
reason there is no competition in this neutral situation is 



because the number of vesicles is infinite. If there were 
a finite number of vesicles then random drift, now act- 
ing at the level of the vesicles, would lead again to the 
fixation of only one vesicle type, in spite of the fact that 
both types have the same survival probability. Use of 
the diffusion approximation, valid in the limit of large N 
and small r such that the product tN is finite, yields a 
simple expression for the fraction of vesicles carrying the 
less fit template A [29j 



X N — 



1 — exp (Nrp) 
1 - exp (Nt) ' 



(5) 



where p is the initial frequency of template A in each vesi- 
cle. Clearly, Yq° = 1 — Y^f. This analytical prediction 
is compared with the results obtained by the numerical 
iteration of the recursions Q in Figure [T] where we have 
used Y° /2 = 1 (hence Y? = if i ^ N/2) so that p = 1/2. 
As expected, there is a good agreement between the ex- 
act numerical and the approximate analytical results for 
small values of the handicap t. In particular, for t = 
we have = p = 1/2. However, we find that 
decreases much faster than exp(— Nt/2) with increasing 

T. 

The point here is to stress that although random drift 
is a key element of models for the evolution of altruistic 
traits, as it enables the fixation of the less fit template in 
some vesicles (provided the handicap r is small), it is a 
serious hindrance to the coexistence of distinct templates 
within a same vesicle. It is in this sense that we can say 
the coexistence issue is more tricky to explain (and actu- 
ally much less studied) than the problem of the altruism 
in nature. 



IV. DETERMINISTIC LIMIT 

Another hindrance to template coexistence, as well as 
to the evolution of altruism via the fixation of the less fit 
template A, is the deterministic pressure in favor of the 
fitter template B that prevails in the limit of infinitely 
large vesicle capacities N — > oo, regardless of the value 
of the cooperation pressure g < 1. In fact, the reason 
that the uniform initial vesicle population (i.e., all vesi- 
cles composed of the same number of templates A and 
B) used to draw Figure [1] resulted in the two antago- 
nistic types of vesicles was the amplification of random 
fluctuations which is an inherent feature of the compet- 
itive dynamics in a finite population. For N — > oo such 
fluctuations are absent and so all vesicles have the same 
composition at all generations. As a result, the compe- 
tition between vesicles is turned off thus leading to the 
fixation of the fitter template in all vesicles (except in the 
degenerate case r = 0, for which the two template types 
coexist forever). 

It is clear then that template coexistence will be possi- 
ble only within a narrow range of the values of the control 
parameters and, in particular, of the vesicle capacity. In 
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FIG. 2: Fraction of vesicles with the two template types in 
their compositions as function of the vesicle capacity for t — 
(i.e., the templates have the same replication rate), g = 0.1 
and m as indicated in the figure. 
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FIG. 3: Critical value of the vesicle capacity as function of 
the ratio between vesicle and template generation times for 
r = and (left to right) g = 0.1, 0.2, . . . , 0.9. Below N c , the 
fixation of one of the template types via random drift bars 
coexistence. 

the following section we determine the regions in the pa- 
rameters space where coexistence is possible. 

V. TEMPLATE COEXISTENCE 

Here we will focus on the fraction fl of vesicles in the 
stationary regime that carry copies of the two template 
types regardless of their redundancies, so that Q = 1 — 

yoo \roo 
1 N 1 ■ 

As pointed out before, in the case of degenerate tem- 
plates r = the sole obstacle to coexistence is the effect 
of random drift, which is intensified for small vesicle ca- 
pacities. This result is illustrated in Figure [21 where we 
show also the effect of increasing the ratio m between the 
vesicle and template generation times. In this, as well as 
in the following figures in which the integer quantities 
N and m are depicted, the continuous lines are simply 
guides to the eye. For N — > oo we find exceptionally that 
— > 1. Two features of this figure are worth emphasiz- 
ing. First, for a fixed value of the cooperation pressure 
g, there is a minimum value of the vesicle capacity below 
which coexistence is unattainable. Second, this threshold 
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FIG. 4: Fraction of vesicles with the two template types in 
their compositions as function of the vesicle capacity for r = 
0.1, g — 0.35 and m as indicated in the figure. 



value, denoted by N c , increases with increasing m, i.e., 
coexistence is inhibited if the local template populations 
evolve faster than the vesicle population. In fact, Figure 
[3] shows that N c increases linearly with m and that the 
slopes of the straight lines decrease in a nontrivial way 
as the cooperation pressure increases. 

We turn now to the analysis of the general case where 
the template types have different replication rates, r > 0. 
In contrast with the degenerate case discussed before, in- 
creasing the capacity of the vesicles N will now favor the 
fixation of the fittest template, thus inhibiting coexis- 
tence. This is exactly what Figured] depicts: the fraction 
of vesicles with the two template types D, is nonzero only 
for a well-defined range of the vesicle capacity, which de- 
creases as m increases. We note that because coexistence 
is impossible in the deterministic regime we have O = 
for N — > oo regardless of the value of m. Figures [5] and 
[5] summarize this finding by showing N c (i.e., the value 
of N at which Q, vanishes) as function of m for several 
values of the control parameters r and g. As hinted in 
Figure |H there is a certain value m = m c (r, g) beyond 
which coexistence is impossible regardless of the vesicle 
capacity. For m < m c , we always find two solutions for 
N c corresponding to the lower and higher vesicle capac- 
ity compatible with coexistence. The lower bound for 
N changes little with variation of g (see Figure [5]) and 
is practically insensitive to variation of t (see Figure [6|, 
provided that the parameter setting remains within the 
coexistence boundary, i.e., m < m c . The upper bound, 
however, is extremely sensitive to variation of those pa- 
rameters. Finally, Figure [7] illustrates how m c depends 
on the parameters t and g. 

Up to now we have emphasized the role of the gener- 
ation times of templates and vesicles as yet another ob- 
stacle to the coexistence of distinct template types, thus 
generalizing previous approaches that assumed that tem- 
plate and vesicle populations were updated (this term is 
appropriate since in both cases it is assumed that gen- 
erations do not overlap) with the same frequency, i.e., 
m = 1 [27l ]. For completeness, in Figure [8] we move m 
to a secondary position and stress the effect of the coop- 
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FIG. 5: Critical value of the vesicle capacity as function of 
the ratio between vesicle and template generation times for 
r = 0.1 and g as indicated in the figure. For fixed m the two 
values of N c determine the range of vesicle capacity within 
which coexistence is possible. 
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FIG. 6: Critical value of the vesicle capacity as function of 
the ratio between vesicle and template generation times for 
g = 0.8 and r as indicated in the figure. 



eration pressure g and handicap r on the coexistence of 
templates A and B. These results illustrate clearly that 
drift is the main obstacle to coexistence in the case the 
handicap r is small, since a large cooperation pressure is 
needed to retain the two template types for small vesicle 
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FIG. 7: Critical value of the ratio between vesicle and tem- 
plate generation times as function of the cooperation pressure 
for r as indicated in the figure. Above m c coexistence is im- 
possible. 
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FIG. 8: Critical value of the handicap parameter as function 
of the cooperation pressure for m = 1 (main graph) and m — 5 
(inset). The vesicle capacities are indicated in the main graph. 
Below t c coexistence is impossible. 



sizes. As r increases, however, the deterministic template 
competition rapidly takes the lead as the main hindrance 
to coexistence. The effect of increasing m, as shown in 
the inset, is to shift non-uniformly the coexistence lines 
to higher values of the cooperation pressure. 



VI. CONCLUSION 

Compartmentalization of unlinked templates is widely 
recognized as a necessary step towards the evolution of 
cellular life, but it is rarely appreciated that confining a 
finite number of templates in a vesicle actually creates a 
new obstacle to coexistence - the fixation of one of the 
template types caused by random drift. Since this dis- 
ruptive effect is enhanced in vesicles of low capacity, one 
should not expect to find such vesicles in a realistic pre- 
biotic scenario. On the other hand, the risk of vesicles of 
large capacity is well-known: the deterministic competi- 
tion between templates results in the fixation of the type 
with the higher replication rate. So very large vesicles 
are not to be expected in a prebiotic scenario too. 

The range of permitted vesicle capacities depends on 
several biologically relevant parameters, and we can spec- 
ulate on the values that produce a sensible scenario. For 
instance, it is now well established that vesicles sponta- 
neously assembled from fatty acid micelles [H], Hf| (see 
also [36|) grow and divide competing for the free micelles 
in the environment. Hence vesicles can do well without 
templates and so the parameter g that appears in the 
survival probability of the vesicles, equation (JTJ, should 
be set to a small value, but not a too small one since 
then it would be impossible to compensate for the pres- 
sures of drift and competition. This observation imme- 
diately excludes low capacity vesicles from our prebiotic 
scenario (see Figure [5J). However, the same figure shows 
that coexistence in high capacity vesicles, say N = 100, 
is possible only if the templates have near degenerate 
replication rates (r <C 1). Of course, one may object 
that this is a precarious situation since sooner or later a 
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mutant template with a higher replication rate will show 
up and destroy this fragile balance. The answer to this 
criticism brings out the main advantage of compartmen- 
talization: the mutant (or parasite) appears in a single 
vesicle which is then unlike to pass the extinction stage 
of the life cycle, since by definition the presence of the 
parasite reduces its chance of survival, equation |T|). The 
parasite is then quickly eliminated as a result of the death 
of the infected vesicle. 

Although one might think it is plausible to assume that 
the templates have a much shorter generation time than 
the vesicles, i.e., m 3> 1, that would prevent coexistence 
altogether (see Figure [7]). In addition, following the arm- 
chair argument given before, one expects the number of 
infected vesicles to increase with the number of template 
replication cycles m simply because mutants appear as 
results of errors during the replication stage. So, though 
counter- intuitive at first sight, we are forced to admit 
that the template generation time must be of the same 
order or even longer than that of the vesicles. Our model, 
as well as all known package models, is not suitable to 
describe this possibility (m < 1), but it seems clear from 
our results that coexistence would be favored in this case. 
Attempt at modifying the present analytical formulation 



to describe the case m < 1 is under way. 

The ultimate goal of theoretical research on prebiotic 
evolution is to come up with a coherent scenario for the 
origin of life. Our study supports the view that such a 
plot begun with a very large population of vesicles capa- 
ble of template-independent reproduction. The acciden- 
tal assimilation of different species of unlinked templates 
that happened to boost the reproduction capability of the 
vesicles assured then that only vesicles containing a spe- 
cial kind of templates - those with near-degenerate repli- 
cation rates and long generation times - would thrive. At 
this point the information crisis was overcome and the 
stage for a genetic takeover was set: the vesicles would 
soon loose their ability to reproduce without the aid of 
templates, and become thus mere vehicles or means for 
template replication. To refine or reject this scenario is 
the main theme of theoretical research on the evolution 
of life [13. 
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